Pulsar Kicks from Active-Sterile Neutrino Transformation in Supernovae 
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Observations of radio pulsars have revealed that they have large velocities which may be greater 
than 1000 km/s. In this work, the efficacy of an active-sterile neutrino transformation mechanism to 
provide these large pulsar kicks is investigated. A phase-space based approach is adopted to follow 
the the transformation of active neutrinos to sterile neutrinos through an MSW-like resonance in the 
protoneutron star to refine an estimate to the magnitude of the pulsar kick that can be generated 
in such an event. Estimates of the magnitude of the pulsar kick created and the overall cooling due 
to sterile neutrino emission suggest that this mechanism can create the large pulsar kicks while not 
overcooling the star. 

PACS numbers: 14.60.St, 97.60.Bw, 97.60.Gb 



I. INTRODUCTION 

Supernova explosions have long been used as a labo- 
ratory for particle physics, especially for the existence 
of light, weakly interacting particles such as axions, ma- 
jorons, and sterile neutrinos. In particular, sterile neu- 
trinos are of interest because one with a mass of several 
keV is a viable - and possibly observable - dark matter 
candidate. The study of the supernova mechanism can 
provide limits on sterile neutrino parameters from the re- 
quirement that the emission of sterile neutrinos from the 
protoneutron star does not overcool the star before the 
emission of active neutrinos which are observed in the 
SN1987A event. 

Astronomical observations show evidence that radio 
pulsars can have very large velocities [ll-[3]. Estimates 
of their average velocities are of the order of a few hun- 
dred km/s, but some pulsars appear to have velocities in 
excess of 1000 km/s 

Radio pulsars are a small subset of the billions of neu- 
tron stars in the galaxy; each neutron star being born in 
a supernova, the cataclysmic end to the life of a massive 
star. It would be natural to look to the supernova explo- 
sion for a mechanism to generate these large observed ve- 
locities of pulsars ( "pulsar kicks" ) . While hydrodynami- 
cal asymmetries in the supernova collapse may be able to 
generate such large velocities, many have failed. Recent 
app roaches are based on two-dimensional simulations 
lll| . but it is unclear whether, in three dimensions, this 
mechanism can account for both the magnitude of the 
kick and the correlation between the spin axis and pul- 
sar velocity that is observed in the data [l2|-[T5|. 

The signature of keV sterile neutrino emission in a su- 
pernova is a significant momentum anisotropy, that gen- 
erates a large pulsar kick ^16^. In addition, the observed 
spin- velocity correlation is a generic result of keV ster- 
ile neutrino emission with a large enough magnitude to 
generate the observed pulsar kicks. The observed pulsar 
kicks can be explained by the emission of sterile neutrinos 
created by active-sterile neutrino transformation within 
the highly magnetized, dense core of a protoneutron star. 



The purpose of this paper is to understand the physics 
of sterile neutrino creation and subsequent emission in 
detail. 

Sterile neutrinos are an intriguing extension of the 
Standard Model toward resolving unexplained phenom- 
ena in astrophysics [l7| . The concept of a sterile neutrino, 
or "right-handed neutrino," as a neutrino that does not 
participate in the weak interaction was introduced by 
Bruno Pontecorvo [ll|. Although the sterile neutrino 
does not directly couple to the strong, weak or electro- 
magnetic interactions, it may couple to the active neutri- 
nos and/or the Higgs boson which may provide a number 
of avenues for sterile neutrinos to be an integral compo- 
nent in astroparticle physics. 

These right-handed neutrinos are natural extensions of 
the Standard Model to account for the fact that neutrinos 
have mass. The seesaw mechanism can allow the mass 
of these right-handed neutrinos to be very large (TsH^H j 
but has also recently been shown that can also allow the 
mass of these neutrinos to be smaller [22| . It is the latter 
of these possibilities that are regarded as sterile neutrinos 
that provide exciting opportunities to explain astrophys- 
ical phenomena. 

This paper will discuss the prospects of a sterile neu- 
trino with a keV mass and vacuum mixing between the 
sterile and active neutrinos. There has been significant 
interest in a keV mass sterile neutrino as a proposed dark 
matter candidate f2^-'28'|, possible mechanism for gener- 
ating pulsar kicks [16, 29, 30], strengthening the super- 
nova explosion |3ll433l |. and as a source of ionization in 
the reionization epoch 3J|. In addition, preliminary as- 
trophysical observations suggest the possibility of the ex- 
istence of a dark matter candidate sterjle neutrino with 
masses of either 5 keV (35i] or 17 keV jSQ] . 

In Section |lT] an active-sterile neutrino transformation 
mechanism to produce pulsar kicks will be presented. 
Section IIIII discusses the calculation of the magnitude of 
the pulsar kick generated, while Section IIVI applies the 
proposed mechanism in a mass element of a protoneutron 
star to estimate the magnitude of the pulsar kick gener- 
ated by active-sterile neutrino transformation. Finally, a 
brief discussion on the mechanism is presented in Section 
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IVl and conclusions are presented in Section |Vj 



II. RESONANT STERILE NEUTRINO 
PRODUCTION 

Let us consider the evolution of neutrino states during 
the first ~ fO s after the formation of a protoneutron 
star. In this discussion, we will consider the simplest 
of active-sterile neutrino mixing between one active and 
one sterile state so that the mass and weak eigenstates 
are related in vacuum by 



= cos 0\Ua) ^ sm ()\iys) 
|j/2) = sin + cos 6*11/5), 



(1) 



where V12 are the mass states, i^a with a — fj,, or 
T is the active neutrino, Vs is the sterile neutrino, and 
9 is the vacuum mixing angle between the active and 
sterile neutrino. As the neutrinos travel through matter, 
the relationship in Eq. ([1]) is slightly altered as the mass 
states become instantaneous mass {i.e., energy) states 
and the vacuum mixing angle is replaced by the effective 
matter mixing angle 9m. The effective matter mixing 
angle can be related to the forward scattering potential 
of Va in matter, Va, by 



sin^ 29„ 



(Jm2)/2p)2sin2 29 



(Jm2/2p)2sin2 29+[{Sm^/2p) cos26'- Va]^' 

(2) 

for a neutrino state with momentum p [25| . For the ster- 
ile neutrino parameters that are relevant in this paper 
(sin^ <C 1, m2 ^ rrii), the difference in the squares of 
the two neutrino masses, Srri^ — m| — m^, can be related 
to the "sterile neutrino mass," rris, by Sm^ sa m^. 

In vacuum, since sin^ 0^1, sterile neutrinos are most 
closely associated with 1^2, the heavier mass state. How- 
ever, as the forward scattering potential increases, so 
does the matter mixing angle. If the forward scatter- 
ing potential increases to a point that is analogous to the 
MSW resonance [stUssI, 



1^ cos 26* 



(3) 



then 9m — 7r/4; the active and sterile neutrino are maxi- 
mally mixed. 

The dynamics of the neutrino state can be discussed 
in terms of three length scales: the neutrino oscillation 
length, the width of the resonance discussed above, and 
the neutrino mean free path. There are two mechanisms 
in the conversion of active neutrinos to sterile neutrinos: 
non-resonant and resonant production of sterile neutri- 
nos. 

Scattering induced, non-resonant production of ster- 
ile neutrinos occurs at a rate proportional to the neu- 
trino scattering rate and sin^ 29m. However, if the neu- 
trino oscillation length is short compared to the scat- 
tering length, then sterile neutrino production will be 
suppressed by the quantum Zeno effect. 



Resonant sterile neutrino production occurs at the 
MSW-like resonances described by Eq. This res- 

onant production mechanism is efficient in converting 
active neutrinos that meet the resonance condition into 
sterile neutrinos as long as the neutrino oscillation length 
at resonance is small compared to the resonance width 
which is small compared to the mean free path. A large 
momentum anisotropy from the conversion of active to 
sterile neutrinos can be produced by this mechanism. 

The neutrino oscillation length is maximized at the 
resonance point. The neutrino oscillation length at reso- 
nance is 



2p 



1? sin 261' 



(4) 



For the adiabatic evolution of these neutrinos, the neu- 
trino oscillation length at resonance should be small com- 
pared to the width of the resonance. 



6t 



1 dVr, 



Va dx 



tan 29, 



(5) 



where d/dx is the spatial derivative along the neutrino 
trajectory. Finally, to describe the coherent evolution of 
the neutrino as it traverses the resonance, the width of 
the resonance should be small compared to the mean free 
path of the active neutrino. 



A 



mfp 



(6) 



where p is the density, mn is the neutron mass, and 
CTweak ~ G\p^ is the weak scattering cross section with 
Gf the Fermi coupling constant [25j . 

If the width of the resonance is small compared to the 
neutrino mean free path, we can approach the evolution 
of neutrinos through the MSW resonances as being coher- 
ent. Neglecting any neutrino magnetic moments, the for- 
ward scattering potential for active neutrinos in a magne- 
tized, neutral and degenerate background can be written 

as [la-Iil 

+ V^{p,B,T-Y,)cos9, (7) 
Vp = -Vo{p)fp{Yi) + V^{p,B,T-Y,)cos9,, (8) 

for f3 — p or T, cos 9^, is the angle between the neutrino 
momentum and magnetic field, where 



fpiYi)^l-Ye- 2Y,„ - 2Y„^ - 2Y,^ - 2Y,^ , (9) 



and 



G o 

yo(p)--|— «3.8pi4 eV, 
V2 riiN 



(10) 



with pi4 the density in units of 10^* g/cm^^. The Yg 
are the local lepton asymmetries and are defined as the 
lepton number divided by the baryon number. 



Yf 



n{ — Tig 
rib 



(11) 
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where rii, and Ub are the number densities of species 
£, its anti-particle and baryons, respectively. For anti- 
neutrinos, the forward scattering potential is the neg- 
ative of the forward scattering potential for neutrinos, 
Va = — Va, and the forward scattering potential for ster- 
ile neutrinos is zero, Vs — 0. 

The terms in the forward scattering potential that de- 
pend on the magnetic field result from forward scattering 
off polarized background electrons and baryons that tend 
to have their spins anti-parallel to the magnetic field. For 



V,^p,B,T;Ye) = i- 



Cn)B, 



(12) 



and for Vn and Vt, 



Vlip,B,T;Ye)^iCe + Cp + Cn)B. (13) 

The difference between the forward scattering potentials 
for Ve and i^^^^r is that Ve have both charged and neu- 
tral current contributions from scattering off electrons, 
while Vp^^r have only neutral current contributions. The 
energy scales in the protoneutron star are well below the 
threshold to thermally populate muon or tau states. 

For degenerate electrons, Ue ^ , so rie « YeUb- The 
neutral current contribution from neutrino scattering off 
polarized electrons is 



,B 



2.0pifBi6 



Ye 

0.04 



1/3 



for scattering off polarized protons is 



CpB ^ 0.85 puBie 



Ye 

0.04 



meV, 



meV, 



and for scattering off polarized neutrons is 
CnB « 15pi4i?i6 meV, 



(14) 



(15) 



(16) 



where Big is the strength of the magnetic field in units 
of 10^6 G El. 

In this paper we will work with 2x2 mixing between 
V/j, and Vs (and their CP-related counterparts). This 
analysis will present a flavor of the physical mechanism 
that drives the pulsar kick through active-sterile neutrino 
transformation. The effects of mixing with all three ac- 
tive species will be discussed in Section HV CI 

At a given time and location in a neutron star, the local 
physical conditions (density, temperature, and magnetic 
field), along with the local chemistry (the Yg), determine 
the resonant neutrino energy as a function of the neutrino 
direction relative to the magnetic field (cos 9^). The reso- 
nance condition, Eq. ([3]), applies to u^^ ^ Vg when > 
or Ufj_ ^ Ps when < 0. 

Sterile neutrinos are created when ^ neutrinos coher- 
ently evolve through a resonance. The resonant momen- 
tum is 



mz 



(17) 



where we have made the assumption that the vacuum 
mixing angle is small, and the tacit understanding that 
Pros > applies to ^ Vs and p^^s < applies to 
^ Vs- It is evident that the resonant momentum 
depends on the local conditions and the direction of the 
neutrino. 

If the neutrino oscillation length at resonance is small 
compared to the width of the resonance, then the coher- 
ent evolution of the neutrino states is adiabatic resulting 
in the efficient conversion of active to sterile neutrinos 
(or anti- neutrinos). That is, all neutrinos traveling in 
a given direction with the resonant momentum will be 
converted; active neutrinos will become sterile, and any 
sterile neutrinos will be converted back to an active neu- 
trino. However, if the neutrino oscillation is of order the 
resonance width, or even larger, the neutrino propaga- 
tion through the resonance is non-adiabatic and as a re- 
sult the active-sterile neutrino conversion is less efficient, 
albeit non-zero. 

v^ or with momenta dictated by the resonance con- 
dition are converted into sterile neutrinos with a proba- 
bility dictated by how adiabatic the conditions are as the 
neutrino states cross the resonance. The width of the 
resonance in momentum space is Picstan20, so a "hole" 
is created in the active neutrino and anti-neutrino dis- 
tribution functions at the resonant momentum with the 
given width. Y^^ changes as neutrinos and anti-neutrinos 
are converted to sterile neutrinos and energy is radiated 
away as the sterile neutrinos freely stream out of the pro- 
toneutron star. 

In the core of the protoneutron star, the degenerate 
neutrons are the dominant component in determining the 
local energy and pressure. The leptons play a smaller 
role in the local energy and pressure budget of which the 
/i-neutrinos and anti-neutrinos have only a fractional in- 
fluence (less so early on when the electrons and electron- 
neutrinos are degenerate while the /i-neutrinos are not). 
So, it would be reasonable to proceed as if the produc- 
tion of sterile neutrinos have no feedback on the physical 
conditions of the neutron star. 

Since the neutrons far outnumber the neutrinos, which 
in turn, far out number the resonant neutrinos that trans- 
form into sterile neutrinos, the concomitant change in 
y^y^is small. The result is that the resonant momentum 
does not change significantly and if the hole in the neu- 
trino distribution function is not repopulated, the sterile 
neutrino production mechanism would stall. However, 
at the high densities found inside the protoneutron star 
core, scattering and neutron bremsstrahlung are efficient 
mechanisms for repopulating the hole created in the neu- 
trino distribution function. 

As neutrinos are created at the resonant momentum, 
they are converted to sterile neutrinos as described above, 
resulting in continued change in YJy^ . Eventually, the 
change in y,y^will become large enough to shift the res- 
onance to thermally populated regions of phase space 
and the process begins anew. The active-sterile neutrino 
transformation resonance sweeps through the active neu- 
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trino distributions converting active neutrinos into sterile 
neutrinos. The sweeping of the resonance continues un- 
til the forward scattering potential is altered so that the 
total active-sterile conversion converts an equal number 
of neutrinos as anti- neutrinos and thus is unchanged. 



resonance width {ires ^ St) and tends toward zero when 
Ires ^ St. Here, we use the Landau-Zener jump proba- 
bility Hi,!!!, 



1 



.-^7/2 



(21) 



III. CALCULATIONS 



which approaches unity for 7 3> 1 and is proportional to 
7 for 7^1. The adiabaticity parameter is 



A. Resonant Production 

As described in the previous section, the chemical evo- 
lution of a mass element in the protoneutron star is dom- 
inated by the scattering of neutrinos into the resonance. 
When z/p are scattered into the resonance and converted 
into Vs, Fi/^ decreases, and when are scattered into 
the resonance and converted into v^, Yiy^ increases. The 
evolution of is important in the determination of the 
spectra of sterile neutrinos emitted as a result of the 
active-sterile neutrino transformation. 

The sterile neutrino production rate is equal to the 
rate active neutrinos are scattered into the hole in the 
neutrino distribution created by resonant active-sterile 
neutrino transformation. The scattering rate of neutrinos 
into the hole in the distribution is proportional to the 
overall scattering rate. In addition, the scattering rate is 
proportional to the size of the resonance, i.e., 



St 



(2^ 



Prcs tan 29 d cos 9^ 



(18) 



where, by symmetry, the azimuthal direction about the 
magnetic field is integrated out. Finally, the scattering 
rate is proportional to the fraction of neutrinos avail- 
able to scatter into the resonant momentum. Thus, the 
sterile neutrino production rate for neutrino states with 
momenta in the direction described by cos 9^, is 



driv, 
"dT 



d cos 9 1, 



Pres tan 29 d cos 9^ 



Pics 



/T2 



47r2 



(19) 

where T^^ is the total scattering rate and V^^^u^ is 
the probability that a with the resonant momentum 
will transform to Vs as it traverses the resonance. 

In the core of the protoneutron star, the neutrino scat- 
tering rate is primarily due to neutrino-nucleon scattering 
[ill, so the total neutrino scattering rate is 



(20) 



At the temperatures of interest, muons are not thermally 
populated, so neutral current interactions are responsible 
for the scattering of the and into the resonant mo- 
mentum states. Consequently, the production rates for 
j/^ — > Us should have the same form as those for Vg 
shown above. 

The probability of conversion should be unity when 
the neutrino oscillation length is small compared to the 



2p 



1 dVu. 



dx 



sin 26*. 



(22) 



The evolution of Y^, is thus dictated by 



dt 



P 



1 \ dt 



dt 



dcos9^. (23) 



The evolution of Y^ is important because the value of 
yjy along with the other lepton numbers and the phys- 
ical conditions in the protoneutron star determine the 
resonant momenta. 

The sterile neutrinos that are created in these reso- 
nances subsequently stream freely out of the protoneu- 
tron star since they no longer interact with the dense en- 
vironment. By symmetry, the total momentum of these 
sterile neutrinos are either parallel or anti-parallel to the 
magnetic field. Thus, in the calculation of the total 
momentum emitted in sterile neutrinos, we need only 
take the component of the neutrino momentum along 
the magnetic field into account when calculating the to- 
tal momentum emitted from the protoneutron star. The 
rate at which momentum is emitted in sterile neutrinos 
per unit mass is 



1 f^^ f dn^^ drijj^ 



dps 

dm dt p J _i \ dt 



dt 



pcos6„ dcos6„. (24) 



The rate at which energy is emitted in sterile neutrinos 
per unit mass is 



dEs 
dm dt 



1 



dn^, 
~~dt 



duj)^ 
"dt 



pdcos 9„. 



(25) 



The calculation presented will consider the sterile neu- 
trino production at a location in a neutron star with a 
fixed value of its density, p and temperature T. The ini- 
tial conditions on the lepton numbers are chosen to be 
Ye = 0.04, Y^^ = 0.07, and F^^, = Y^^ = 0. 

With the initial conditions listed above, is positive, 
which means that is negative for all allowed values of 

cos 9^ as long as B < 8pi{'^ x 10^^ G. The form of the 
resonance condition, Eq. (|17p . tends to have higher reso- 
nant momenta for neutrinos directed into the hemisphere 
parallel to the magnetic field (hereafter, the "northern 
hemisphere") than the hemisphere anti-parallel to the 
magnetic field (hereafter, the "southern hemisphere"), 
resulting in an anisotropy of the emitted momentum in 
sterile neutrinos. 
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B. Reabsorption 

The same MSW-like resonance that transforms active 
neutrinos into sterile neutrinos will also transform sterile 
neutrinos back into active neutrinos if the neutrino's mo- 
mentum is once again resonant at a later point in time 
along its worldline. As a consequence, the total momen- 
tum in sterile neutrinos that escape from the neutron 
star is reduced as sterile neutrinos that would otherwise 
escape the star are converted back into active neutrinos 
and are rethermalized in the sea of active neutrinos in 
another location. The result is that evolution of Y^^ is 
altered in a manner that depends on the location of the 
mass element in the star and on the regions of phase 
space in the v^j_ and/or distribution functions that are 
reabsorbed. 

In a mass element of the star, the sterile neutrinos that 
are emitted in the hemisphere directed out of the star will 
not traverse another resonance and will thus escape from 
the star. On the other hand, the sterile neutrinos that 
are emitted in the hemisphere directed toward the center 
of the star will encounter another resonance and thus be 
reabsorbed in another mass element on the other side of 
the star. The many-decades range of densities (and rela- 
tively small range of ) in the core of the protoneutron 
star effectively assures that inward emitted neutrinos will 
encounter another resonance. (More general geometries 
may result in multiple resonances along the trajectory of 
a sterile neutrino. In such a case, the interior resonances 
will have a reabsorption rate comparable to the emission 
rate because there will be resonant neutrinos in both in- 
ward and outward directions. Thus, only the outermost 
resonances will contribute to the overall effect.) 

Roughly, the global kinematic effect of this descrip- 
tion is that only the sterile neutrinos emitted out of the 
star escape from the star because these neutrinos are 
not reabsorbed. On the other hand, the local effects on 
energy transport and the evolution of Yj/^^ depend only 
on the sterile neutrinos emitted into the star, because 
a similar number of neutrinos are emitted in an out- 
ward direction as are reabsorbed as neutrinos traveling 
in the same direction, but have traversed the center of 
the star. To properly determine the transport of energy 
from one location in the protoneutron star to another 
through the transformation, i/^ — > ;/s ^ j/^, is a complex, 
multi-dimensional ray tracing transport problem. (Refs. 
[3^ [33I considered a similar energy transport scheme in 
one dimension.) 

In a mass element of the protoneutron star, the number 
of sterile neutrinos that are reabsorbed is roughly half the 
number of sterile neutrinos emitted. Thus, it is a good 
approximation to halve the rate of change of Y^,^ in Eq. 
(|23p . The momentum and energy rates emitted in sterile 
neutrinos can be estimated by averaging the output of 
diametrically opposed mass elements. Consequently, the 
momentum and energy emission rates in Eqs. (j24p and 
are also roughly halved. The results in the following 
section will be presented with respect to these reabsorp- 



tion effects on the chemical and kinematic evolution and 
will employ the approximations stated above. 



IV. RESULTS 
A. An Example 

As an example of the the sterile neutrino production 
mechanism discussed earlier, we focus on a mass element 
in a neutron star with T = 20 MeV, and B = 10^'^ G. 
The initial conditions stated in the previous section on 
the lepton numbers are also employed: Ye = 0.04, 
Y^^ = 0.07, and Y^^^ = Y^^ = 0. Finally, the density 
scale height, \d\ogV^/ dr\^^ ~ 10 km, is used in calculat- 
ing the adiabaticity parameter. 

Throughout this section, we consider the parameters 
of the mass element in the protoneutron star remain un- 
changed throughout the calculation. While this may be 
a good assumption on the coUisional timescale, as men- 
tioned in the previous section, on the longer dynamical 
timescales considered in the calculation, these assump- 
tions ought to break down. A discussion is presented in 
Section IIV CI to qualitatively discuss consequences of a 
multitude of effects that may be present if this mecha- 
nism were to be included in a self-consistent simulation 
of the formation of a pulsar. 

Figure [T] summarizes the chemical and kinematic evo- 
lution due to the production of sterile neutrinos of a 
mass element in the protoneutron star with these phys- 
ical characteristics and a density of 2 x 10^^ g/av?. 
The upper left graphs shows the evolution of Y^^ . The 
two graphs on the right show the rate per unit mass 
of the momentum of the emitted sterile neutrinos (up- 
per) and energy of sterile neutrinos emitted (lower). 
The lower left graph shows the ratio of these momen- 
tum and energy rates, a measure of the asymmetry in 
the sterile neutrino emission. If this ratio [dps/dEs — 
{dps / dm/ dt) / {dEs / dm/ dt)] is equal to -1, then the ster- 
ile neutrinos are emitted solely opposite to the direction 
of the magnetic field. 

The general shapes of the curves in Figured] are conse- 
quences of the mathematical form of the forward scatter- 
ing potential and the resonances in this problem. There 
are three generic features of the sterile neutrino produc- 
tion mechanism discussed in this paper. (1) Y^, increases 
from its initial value of and asymptotes to a fixed value 
that depends on the initial value of the lepton numbers in 
electrons and electron neutrinos. (2) The net momentum 
of the emitted sterile neutrinos is opposite the direction of 
the magnetic field. (3) The magnitude of the momentum 
and energy emission rates is initially large, but is sup- 
pressed by multiple orders of magnitude on a timescale 
of a fraction of a second. 

These characteristics of the active-sterile transforma- 
tion mechanism in the protoneutron star core can be 
explained by understanding the regions in phase space 
where the transformation is taking place. Figure [5] 
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FIG. 1: Results for the calculation of the production of sterile neutrinos in a mass element of the protoneutron star. Upper 
left panel is the temporal evolution of Y,^^. Upper right panel is the rate per unit mass of momentum of sterile neutrino 
produced in the opposite direction of the magnetic field. Lower left panel is the ratio of the rates of net momentum along the 
magnetic field to energy emitted in sterile neutrinos. Lower right panel is the rate per unit mass of energy emitted in sterile 
neutrinos. Sterile neutrino parameters: nis — 5 keV, sin^ 29 — 10~^. Physical parameters: Big = 10, T = 20 MeV, pi4 = 0.02, 
\d\ogV^/dr\''^ = 10 km. Initial conditions: Ye = 0.04, Y^^ = 0.07, Y^^ = Y^^ = 0. 



shows the evolution of three important resonant mo- 
menta through the first half-second of the process, p+i 
(p-i) represents the resonant momentum for neutrinos 
that are traveling parallel (anti-parallel) to the magnetic 
field. (Thus, cos9^ = -|-1(— 1).) In a similar vein, po 
represents the resonant momentum of neutrinos that are 
traveling perpendicular to the magnetic field. The lower 
figure is the sterile neutrino production rate in a logarith- 
mic scale with each tic mark representing a difference in 
the rate of two orders of magnitude. The basis of the 
emission asymmetry can be seen as the sterile neutrino 
rate for neutrinos emitted into the southern hemisphere 
is orders of magnitude larger than the rate for those emit- 
ted into the northern hemisphere. 

In Figure [21 i^^ are converted to Ds into the south- 
ern hemisphere with momenta between p_i and pq. On 
the other hand, 9^ are converted to Ps into the northern 
hemisphere with momenta between po and p+i- 

increases because are transformed into sterile 
neutrinos, while the z/^ are not affected by this resonant 
conversion. As Yi,^ increases, f^{Yi) decreases, resulting 
in an increase in each of the respective resonant mo- 
menta. As the resonant momenta shift higher, the sterile 
neutrino production rates become suppressed both expo- 
nentially due to the paucity of active neutrinos with these 



high energies available to scatter into the resonance and 
linearly as the conversion becomes non-adiabatic. 

The asymmetry in the overall momentum of the emit- 
ted sterile neutrinos is because p-i < p+i. As is demon- 
strated in Figure [21 the lower resonant momenta cor- 
respond to much larger sterile neutrino emission rates. 
Thus, there are more sterile neutrinos emitted into the 
southern hemisphere than the northern hemisphere, re- 
sulting in a net momentum pointing anti-parallel to the 
magnetic field. As are converted into sterile neutri- 
nos, causing the increase in Y^^as discussed previously, 
the resonant momenta become larger where the produc- 
tion rate decreases exponentially resulting in the steep 
decline in the rates of momentum and energy emitted in 
sterile neutrinos as seen in Figure [1] 

B. Generating the pulsar kick and cooling the star 

Figure [3] illustrates the energy emission per unit mass 
of sterile neutrinos as a function of density assuming 
T = 20 MeV and Sie = 10. The curves are the total en- 
ergy emitted in the first half-second, the first second, and 
the first three seconds. It can be seen that most of the 
sterile neutrino emission occurs early on in the process. 
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FIG. 2: The evolution of resonant momenta p+i, p-i and 
po (upper graph) using the same parameters as in Figure [1] 
p+i (p-i) is the resonant momentum for neutrinos traveling 
parallel (anti-parallel) to the magnetic field, po is the reso- 
nant momentum for neutrinos traveling perpendicular to the 
magnetic field. The lower graph is the sterile neutrino pro- 
duction rate (each tic represents two decades difference in the 
production rate). 
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FIG. 3: The energy per unit mass of sterile neutrino emission 
as a function of density integrated over the first 0.5 s (solid), 
1 s (dashed), and t = 3 s (dot-dashed). All the other physical 
and neutrino properties and initial conditions are the same as 
in Figure [1] 



P14 

FIG. 4: The total momentum per unit mass of sterile neu- 
trinos emitted from a mass element in the protoneutron star 
as a function of density at times t = 0.5 s (solid), t = 1 s 
(dashed), and t = 3 s (dot-dashed). All the other physical 
and neutrino properties and initial conditions are the same as 
in Figure [T] 



is ^ 6 X 10^^ erg and the resultant pulsar kick velocity 
is ~ 2000 km/s, opposite the direction of the magnetic 
field. 

While the strength of the pulsar kick is larger than 
necessary to explain the observed pulsar velocities, there 
are other effects - discussed below - that may decrease 
the magnitude of the pulsar kick. The observation of 
active neutrinos from SN1987A implies that the total en- 
ergy emitted in sterile neutrinos should be ^ 10^^ erg 
|48l |49| , lest the emission of sterile neutrinos overcool 
the protoneutron star, reducing the energy budget left 
for active neutrino emission to match observations. It 
appears that the active-sterile neutrino transformation 
mechanism presented can generate a sufficient kick to 
explain the observational data, yet not overcool the pro- 
toneutron star. 



C. Other Considerations 



as longer integration times result in disproportionately 
smaller increases in the overall energy emission. 

Figure 13] shows the output in terms of the net momen- 
tum in sterile neutrinos emitted from the protoneutron 
star. The curves show the total momentum per unit mass 
that is emitted in the form of sterile neutrinos from the 
protoneutron star as a function of the local density at 
three different times in the evolution. Once again, it can 
be seen that most of the sterile neutrino emission occurs 
early on in the process. The momentum output is nega- 
tive because the net momentum emitted is anti-parallel 
to the magnetic field as discussed above. 

To determine the overall magnitude of the effect, the 
results from Figures |3] and U] must be integrated over a 
density profile for the protoneutron star. Estimating the 
density profile from Ref. [l^l, the total energy emitted 



The calculation performed above involved the determi- 
nation of the net momentum emitted in sterile neutrinos 
from a volume element in a protoneutron star. Although 
the reabsorption of sterile neutrinos due to subsequent 
traversal of resonances is a complex process, the calcula- 
tion presented above attempted to account for this effect 
by self-consistently reducing the sterile neutrino emission 
rates. However, there are three other considerations to 
take into account when determining the magnitude of 
the active-sterile neutrino transformation mechanism for 
large pulsar kicks: dynamical evolution of the physical 
properties of the mass element, the dynamical evolution 
of the active neutrino seas, and sterile neutrino mixing 
with all three active fiavors. 

To accurately model the momentum generated in 
sterile neutrinos, this mechanism for generating sterile 
neutrinos through active-sterile neutrino transformation 
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would need to be calculated alongside a supernova sim- 
ulation. The estimates presented in this paper suggest 
that this mechanism can provide a sufficient kick yet not 
overcool the star. 

The interplay between the emission of a significant 
fraction of the binding energy of the star and the dy- 
namical evolution of the collapsing star can, in turn, af- 
fect the total sterile neutrino output. The temperature 
of the mass element has a strong leverage on the output 
of sterile neutrinos. By inspection of the sterile neutrino 
production rate, the total energy output rate in sterile 
neutrinos is proportional to T^-e ^5 jjjost of the ster- 
ile neutrinos are produced non-adiabatically and 6 if they 
are produced adiabatically). Thus, small changes in tem- 
perature can result in significant changes in the energy 
and momentum emission rates. However, one must be 
cognizant not to overcool the star with respect to the in- 
ferred temperature of the neutrinosphere from SN1987A 

Even without the emission of sterile neutrinos, the tem- 
perature of a mass element in the protoneutron star will 
evolve. At the densities of interest for the calculation pre- 
sented above (just within the neutrinosphere), the tem- 
perature will decrease on a roughly 10 second timescale 
as the heat generated by the infall and bounce stages 
of the supernova are both emitted through the neutri- 
nosphere and diffuse toward the center of the protoneu- 
tron star [5^ . The emission of sterile neutrinos will affect 
the temperature profile of the protoneutron star at late 
times as the total energy emitted becomes a more sig- 
nificant fraction of the total energy budget. In addition, 
the mechanism that generates the large magnetic fields 
necessary for this effect could decay as energy is removed 
by sterile neutrino emission. 

Most of the sterile neutrino emission occurs in the ffist 
^ 50 ms. In half a second, roughly 80% of the total mo- 
mentum is emitted (assuming unchanging temperature 
and magnetic field). The favorable conditions of temper- 
ature and magnetic field presented in this work need only 
be present for a fraction of a second in order to produce 
a kick consistent with observations. Hence, the evolution 
of the temperature and magnetic field strength is unlikely 
to affect the strength of the pulsar kick generated by this 
mechanism beyond the tens of percent level. 

A number of dynamical effects can alter the results 
presented here for a static region. As the protoneutron 
star collapses, the density scale height throughout the 
star will be reduced. As a result, the adiabaticity pa- 
rameter will proportionally shrink, reducing the number 
of sterile neutrinos created, especially at the highest ener- 
gies. The dynamical timescale is likely to be similar to or 
longer than the timescale for changes in the temperature 
and magnetic field, so these dynamic effects will reduce 
the strength of the kick, albeit less so than changes in 
the temperature and magnetic field. Gradients in den- 
sity and temperature will result in gradients in Y^, which 
will lead to diffusion that will act to reduce the disparity 
of 5^1/, in the protoneutron star. However, the neutrino 



diffusion timescale is likely longer than the peak sterile 
neutrino emission timescale, so the effects of this diffu- 
sion are minimal in changing the mechanism. 

In this calculation, it was assumed that there was 
only vacuum mixing between the sterile neutrino and 
j/p. However, if the sterile neutrino mixes with one ac- 
tive neutrino, it would likely also mix with the other ac- 
tive neutrinos as well. Each active neutrino species will 
undergo active-sterile neutrino transformation in a sim- 
ilar fashion to the mechanism described for i/^. Each 
will exhibit a similar asymmetry in the momentum of 
emitted neutrinos resulting from the term in the forward 
scattering potential that depends on both the strength 
and direction of the magnetic field. All three flavors 
of active anti-neutrinos will be resonantly transformed 
into sterile neutrinos leading resulting in the increase in 
each respective lepton number. Since each flavor's lepton 
number will increase, the associated increase of the res- 
onant momenta (as illustrated in Figure [2]) will be more 
rapid. The result is the combinations of two differing 
effects: strengthened by all three flavors being involved 
in the transformation, yet weakened because the sterile 
neutrino emission rates are suppressed as the resonant 
momenta increase. 

If j/g mixes with the sterile neutrino, then a delayed 
kick, urea production mechanism can occur (29l. [30|. In 
this urea production mechanism, the standard urea re- 
actions, + p ^ n + i/g and e'^ + n p + 9^, can 
produce sterile neutrinos instead of electron neutrinos 
with a rate suppressed by sin^20„j. The neutrino states 
produced will have an asymmetry with more neutrinos 
created anti-parallel to the magnetic field than parallel 
to the field [5l|. Thus, any contribution from the urea 
mechanism will augment a pulsar kick generated by the 
active-sterile neutrino transformation mechanism. 

In the urea mechanism, non-resonant active-sterile 
neutrino transformation in the core requires a few sec- 
onds for V"e — >■ 0. Once the forward scattering poten- 
tial for Ve reaches zero, the non-resonant sterile neu- 
trino production rate is maximized f29', 's^ . The mag- 
nitude of this effect depends the conditions in through- 
out the core of the protoneutron star many seconds af- 
ter the active-sterile neutrino transformation mechanism 
was initiated. The large emission of energy in sterile 
neutrinos during the active-sterile neutrino transforma- 
tion mechanism phase will affect the temperature and 
magnetic field strength throughout the core, reducing 
the available energy budget for further emission in sterile 
neutrinos. 



V. DISCUSSION 

The active-sterile neutrino transformation mechanism 
discussed in this paper can feasibly provide the requi- 
site pulsar kick to explain the highest observed pulsar 
velocities. It is interesting that the phase space of sterile 
neutrino parameters (mass and mixing angle with active 
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neutrinos) where active-sterile neutrino transformation 
can explain pulsar kicks overlaps with the phase space 
of sterile neutrino parameters that provide a viable dark 
matter candidate (see Ref. [izl)- It is enticing that two 
astrophysical phenomena could be explained by the same 
sterile neutrino, a straight-forward extension to the Stan- 
dard Model. Combined with tantalizing preliminary ev- 
idence for sterile neutrino dark matter that is consistent 
with this overlap in phase space [35| makes for an inter- 
esting window in sterile neutrino physics. 

Although several mechanisms propose to explain pul- 
sar kicks without the introduction of sterile neutrinos, 
there are deficiencies that plague these mechanisms. The 
evolution of close binaries [52l | and the asymmetric emis- 
sion of electromagnetic radiation jssj have been proposed 
to solve the pulsar kick problem, but neither provide the 
requisite kick to explain the largest pulsar velocities. 

Hydrodynamic mechanisms for generating pulsar kicks 
during the supernova collapse can result in pulsar ve- 
locities of a few hundred km/s, comparable to the av- 
erage observed velocities, but it remains to be seen if 
they can explain the largest observed pulsar velocities 
[sl-fllj. In addition, it is unclear whether the observed 
correlation between pulsar velocity and spin axis [T2!- 
[Tsj can be explained with hydrodynamic mechanisms in 
three-dimensional simulations. The active-sterile neu- 
trino transformation mechanism discussed in this paper 
would be consistent with the velocity-spin correlation be- 
cause although the generated pulsar kick is parallel to 
the magnetic field, any discrepancy between the mag- 
netic field and protoneutron star spin will be averaged 
out due to the rotation of the star . 

The considerations presented in this work do not 
change the general range in sterile neutrino parameter 
space (mass and mixing angle) from past work on this 
mechanism. Requiring that the evolution of the neutri- 
nos are both coherent and adiabatic constrain the ster- 
ile neutrino parameter space available to provide pulsar 
kicks through resonant active-sterile neutrino transfor- 
mation. This parameter space can be found in Refs. 

This work can put some basic constraints on the phys- 
ical conditions within the core of the protoneutron star 
where the active-sterile neutrino transformation is taking 
place (10^^ < pi4 ^ 1). Figure[S]details the magnitude of 
the pulsar kick and the energy emitted in sterile neutrinos 
as a function of the magnetic field-temperature parame- 
ter space for these densities of interest in the protoneu- 
tron star. In order for the active-sterile neutrino trans- 
formation mechanism to produce the pulsar kicks, for a 
given magnetic field strength, the temperature must be 
large enough to produce the fastest kicks 1000 km/s) 
yet small enough that the total energy emission in sterile 
neutrinos does not overcool the star (< 10^'^ erg). 

The active-sterile neutrino transformation mechanism 
for generating pulsar kicks requires large magnetic fields, 
~ 10^^~^^ G, to generate the largest observed pulsar 
kicks. These requisite magnetic fields need only be 
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FIG. 5: Estimates of the total energy emitted in sterile neu- 
trinos and the magnitude of the pulsar kick created by active- 
sterile neutrino transformation mechanism in the parameter 
space of magnetic field (-Bie) and temperature. The pro- 
toneutron star parameters are those in the relevant densities 
(10~^ < pi4 < 1) for the mechanism. The solid lines are con- 
tours of the magnitude of the pulsar kick, while dashed lines 
are contours of the energy emission. 



present for a fraction of a second in the lifetime of the col- 
lapsing protoneutron star for the mechanism to provide a 
sufficient pulsar kick. The magnetic fields may be large, 
as the virial theorem allows for magnetic field strengths of 
~ 10^® G or more [s^ - lsoj . Other considerations suggest 
that it is plausible that these large magnetic fields may 
occur inside some protoneutron stars for the timescales 
needed to generate the pulsar kick as described in this 
work 15614591. 



VI. CONCLUSIONS 

Ref. [m introduced the notion that if there was a ster- 
ile neutrino with a mass of order a few keV with a vacuum 
mixing with active neutrinos that pulsar kicks could be 
explained by resonant active-sterile neutrino transforma- 
tions in a supernova explosion. The authors argue that 
anisotropic emission of sterile neutrinos can create an 
asymmetry in the momentum emitted in sterile neutrinos 
which is a few percent of the total momentum emitted in 
sterile neutrinos. The focus of this paper was to elabo- 
rate upon this mechanism toward providing pulsar kicks 
that can explain the large velocities of pulsars that have 
been observed. 

In this work, a phase space-based approach was taken 
to describe the active-sterile neutrino transformation. At 
any point in the star, a shell in phase space of active 
neutrinos have an MSW-like resonance and will trans- 
form into sterile neutrinos with an efficiency that is dic- 
tated by the adiabaticity of the transition. If the width 
of the resonance is large compared to the neutrino os- 
cillation length, the transition is adiabatic and active 
neutrinos are efficiently converted into sterile neutrinos. 
On the other hand, as the width of the resonance gets 
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smaller, the transition is less adiabatic and the trans- 
formation efficiency suffers. The details of this transfor- 
mation process depend on the density and temperature 
profiles within the neutron star and the magnitude and 
direction of the magnetic field. 

The result is that when the resonant conversion of 
active neutrinos to sterile neutrinos occur at a loca- 
tion within the star where the scattering rate is high, 
a large flux of sterile neutrinos can be created. The 
form of the forward scattering potential, and thus the 
resonant momenta, creates an asymmetry in the emitted 
momentum of sterile neutrinos generated by the resonant 
active-sterile neutrino transformation. For realistic val- 
ues of density, temperature and magnetic field strength, 
an asymmetry on the order of tens of percent can be 
exhibited in the emitted sterile neutrinos. With such a 
large asymmetry, enough sterile neutrinos may be emit- 
ted from the core of the supernova to provide the pulsar 
kick, yet not so many that the sterile neutrinos will over- 
cool the star, contradicting the active neutrino observa- 
tions from Supernova 1987A. 

However, there are a number of other considerations 
that are beyond the scope of this work, but are impor- 
tant to determine the magnitude of the pulsar kick and 
associated cooling by sterile neutrinos. While the esti- 
mates presented in this work provide hopeful suggestions 



that the active-sterile neutrino transformation mecha- 
nism can produce the requisite pulsar kick without over- 
cooling the star, it is necessary for these calculations to 
be made alongside supernova simulations to determine 
the strength of the pulsar kick and the magnitude of the 
cooling due to the emission of sterile neutrinos. 

While the existence of a light sterile neutrino may con- 
tinue to be an unresolved problem in particle physics, 
the possibility of a light sterile neutrino presents exciting 
opportunity in astrophysics to resolve unexplained phe- 
nomena such as dark matter and pulsar kicks. Although 
further exploration is necessary to determine the magni- 
tude of the effect, active-sterile neutrino transformation 
in a supernova presents an intriguing possible mechanism 
toward the explanation of such large velocities observed 
in pulsars. 
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